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Abstract 

We quantify the degree of nonlinearity and stochasticity of the clustering of biased objects, using 
cosmological N-body simulations. Adopting the peaks and the halos as representative biasing models, we 
focus on the two-point correlation of the biased objects, dark matter and their cross-correlation. Especially, 
we take account of the effect of redshift-space distortion and attempt to clarify the scale-dependence and 
the time-dependence by analyzing the biasing factor and the cross-correlation factor. On small scales, 
stochasticity and nonlinearity become appreciable and strongly object-dependent, especially in redshift 
space due to the pair-wise velocity dispersion of the biased objects. Nevertheless, an approximation of 
deterministic linear biasing 5 abj ~ & ob j5 ma3B works reasonably well even in the quasi-linear regime r lOh^ 1 
Mpc, and linear redshift-space distortion explains the clustering amplitudes in redshift space in this regime. 
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Introduction 

Luminous objects, such as galaxies and quasars, are not direct tracers of the mass in the universe. In fact, the 
; I ■ difference of the spatial distribution between luminous objects and dark matter, or the bias, has been indicated 
| from a variety of observations. In order to confront theoretical model predictions for the mass distribution against 
observational data, one needs a relation of density fields of mass and luminous objects. Consider the density contrasts 
of visible objects and mass, 5 ohj (x,z\R) and <5 maB3 (x, z\R), at a position x and a redshift z smoothed over a scale 
R. In general, the former should depend on various other auxiliary variables A defined at different locations x' 
and redshifts z smoothed over different scales R' in addition to the mass density contrast at the same position, 
{ m „ s (x, z\R). While this relation can be schematically expressed as 

5 obJ (x, z\R) = ^[x, z, R, <S maaB (x, z\R),A{x.', z'\R'), . . .}, (1) 

it is impossible even to specify the list of the astrophysical variables A, and thus hopeless to predict the functional 
form in a rigorous manner. Therefore if one simply focuses on the relation between <5 obj (x, z\R) and <5 ma3B (x, z\R), the 
relation becomes inevitably stochastic and nonlinear due to the dependence on unspecified auxiliary variables A. 

For illustrative purposes, define the biasing factor as the ratio of the density contrasts of luminous objects and 
mass: 

(S obj (x,z|fl) z R, <5 ma3B (x, z\R), z'\R'), . . .] 

B *><**W = 5ZMR) = S^(x,z\R) • (2) 

Only in very idealized situations, the above nonlocal stochastic nonlinear factor in terms of <5 ma3B may be approximated 

by 
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1. a local stochastic nonlinear bias: 



B obj (x, z\R) = 6i s b °' [x, z, R, <5 ffla88 (x, z|i?), X(x, z\R), ...], 



(3) 



2. a ZocaZ deterministic nonlinear bias: 



B obi (x, = 6<£>[z, fl, <5 maaB (x, z\R% 



(4) 



and 



3. a local deterministic linear bias: 



B obJ (x,z|i?) = & obJ (z,i?) 



(5) 



From the above point of view, the local deterministic linear bias is obviously unrealistic, but is still a widely 
used conventional model for biasing. In fact, the time- and scale-dependence of the linear bias factor &obj(z, R) was 
neglected in many previous studies of biased galaxy formation until very recently. Currently, however, various models 
beyond the deterministic linear biasing have been seriously considered with particular emphasis on the nonlinear 
and stochastic aspects of the biasing (Pen 1998; Tegmark, Peebles 1998; Dekel, Lahav 1999; Taruya, Koyama, Soda 
1999: Tegmark, Bromley 1999; Blanton et al. 1999; Somerville et al. 1999; Blanton et al. 2000; Taruya 2000; Taruya, 
Suto 2000; Yoshikawa et al. 2000). 

In the present paper, we examine two fairly simplified but well-defined popular models for biasing, i.e., peaks 
and halos, paying particular attention to their stochastic and nonlinear nature on the basis of cosmological N- 
body simulations. Especially we take account of the redshift-space distortion effect on biasing which has not been 
studied in previous work. Our main purpose is to find the validity and limitation of the local deterministic linear 
approximation 5 ohi — b abj S mMS for the general stochastic and nonlinear biasing. By analyzing the two-point statistics 
of the biased objects and the mass distribution, we quantitatively measure the stochasticity and nonlinearity and 
investigate the scale-dependence and the time-evolution of clustering amplitude both in real and redshift spaces. 

The nonlinear stochastic biasing for two-point statistics is not a simple theoretical issue, but important in the 
proper comparison with observation. Although the biasing properties for one-point statistics have been extensively 
discussed in literature using numerical simulations and/or analytic models, a similar investigation of biasing in 
two-point statistics has not yet performed. This is because nonlinear stochastic biasing was originally formulated 
using the one-point statistics. In this paper, we quantify the biasing features of clustering in terms of the two-point 
correlation function. In particular, we focus on the biasing factor 6 obj and the cross correlation coefficient r obj defined 
by equations. (0) and (||). Those quantities are related and compared to the biasing parameters defined in terms of 
the one-point statistics. 

The paper is organized as follows. In section 2, we briefly describe the simulated catalogues and evaluate the 
two-point correlation functions for the peaks and dark halos. We then examine the scale-dependence and the time- 
evolution of biasing measured from those statistics in section 3. In particular we pay attention to the difference of 
the biasing properties of peaks and halos. The results are compared with simple theoretical models and predictions 
from linear theory of redshift-space distortion. In addition, we calculate the variances of peaks, halos and dark 
matter particles taking into account the point process. Section 4 is devoted to the conclusions and discussion. 

2. Peak and halo catalogues from cosmological N-body simulations 

2.1. Simulations 

In this paper, we use a series of cosmological N-body simulations in cold dark matter (CDM) cosmogonies (Jing, 
Suto 1998), whose parameters are listed in Table 1. All the models employ N = 256 3 dark matter particles in the 
periodic comoving cube of the boxsize Lbox = 300/i _1 Mpc, and are evolved on the basis of the Particle-Particle - 
Particle-Mesh (P 3 M) method. The initial conditions of the particle distribution match the CDM transfer function of 
Bardeen et al. (1986) characterized by the shape parameter, T. The RMS mass fluctuation amplitude at 8/i -1 Mpc, 
og, is normalized according to the cluster abundance (Kitayama, Suto 1997). We analyze three realizations for each 
cosmological model at z = and z — 2.2, and one realization for LCDM model at z = 0.6, 1.0, and z = 3.4 to 
examine the time evolution. 

As specific models of biasing, we consider two popular cases; primordial density peaks and dark matter halos. 
The density peaks are selected according to the algorithm of Mo, Jing & White (1997). The initial density field 
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is smoothed with a Gaussian window. The filtering radius of R[ — 0.54 h~ x Mpc is adopted for the smoothing so 
that the typical size of the peaks becomes relevant for galactic-sized objects (Bardeen et al. 1986). Given a density 
field and a threshold peak height, this algorithm predicts the number of the peaks that each simulation dark matter 
particle will carry. The number of peaks per dark matter particle in our simulations is always less than 1. So we 
have picked up the peaks by randomly selecting simulation particles with the selection probability set equal to its 
number of the peaks for each particle. This way guarantees that a correct total number of peaks is reproduced for 
the Gaussian random field. Those peak particles subsequently move according to the gravitational field, and thus 
their total number is conserved. We choose ^th = 1.0, 2.0, and 3.0 as the threshold of the peak height v (Table 2). 
As for the dark matter halos, these are identified using the standard friend-of-friend algorithm with a linking length 
of 0.2 in units of the mean particle separation. We select halos of mass larger than the threshold Mth = 2.0 x 10 12 , 
5.0 x 10 12 and 1.0 x 10 13 hr x M Q (Table 3). The halos having the dark matter particle iV mass < 10 have been 
excluded. 

Figure 1 depicts the distribution of dark matter particles (upper-panel), peaks (middle-panel) and halos (lower- 
panel) in LCDM model at z = and z = 2.2 within a circular slice (comoving radius of 150/i _1 Mpc and thickness 
of 15ft. _1 Mpc). We locate a fiducial observer in the center of the circle. Then the comoving position vector r for a 
particle with a comoving peculiar velocity v at a redshift z is observed at the position s in redshift space: 



where H (z) is the Hubble parameter at z. The right panels in Figure 1 plot the observed distribution in redshift space, 
where the redshift-space distortion is quite visible; the coherent velocity field enhances the structure perpendicular 
to the line-of-sight of the observer (squashing) while the virializcd clump becomes elongated along the line-of-sight 
(finger- of- God). Although the distribution of peaks in Figure la and lb seems rather different, this is simply due 
to the strong concentration along the filamentary regions in Figure la, and the total number of peaks is the same. 
Figure 1 qualitatively illustrates that the redshift-space distortion is sensitively dependent on the objects, and this 
is what we will attempt to quantify statistically in the rest of paper. 

2.2. Two-point statistics 

We use the two-point correlation functions to quantify stochasticity and nonlinearity in biasing of peaks and halos, 
and explore the signature of the redshift-space distortion. Since we are interested in the relation of the biased objects 
and the dark matter, we introduce three correlation functions; the auto-correlation functions of dark matter and the 
objects, £ mm and £ OQ , and their cross-correlation function £ om . In the present case, the subscript o refers to either 
h (halos) or v (peaks). We also use the superscripts R and S to distinguish quantities defined in real and redshift 
spaces, respectively. We estimate those correlation functions using the standard pair-count method. The correlation 
function is evaluated under the distant-observer approximation. 

Those correlation functions are plotted in Figures 2 and 3 for peaks and halos, respectively. Their qualitative 
behavior in redshift space is already well-known and easy to understand (e.g, Kaiser 1987; Hamilton 1998). The 
correlation functions of biased objects generally are larger than those of mass amplitudes. In nonlinear regimes 
(£ > 1) the fmger-of-god effect suppresses the amplitude of £( s ) relative to £*- R \ while £( s ) is larger than £( R ) in 
linear regimes (£ < 1) due to the coherent velocity field. In the next section, we present more quantitative analysis 
in terms of the stochastic nonlinear biasing scheme and compare with the existing theoretical predictions valid in 
linear regimes. 

3. Scale-dependence and time-evolution of biasing parameters 

In order to quantify the behavior of biasing, we introduce the biasing factor and the cross-correlation factor: 




1 



(8) 



(7) 



,mm 



following their counterparts defined in terms of the one-point statistics (e.g., Pen 1998; Tegmark, Peebles 1998; Dekel, 
Lahav 1999; Taruya, Koyama, Soda 1999; Taruya 2000). In equations (0) and (||), 6 obj measures the relative strength 
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of the clustering of objects and r obj characterizes the degree of nonlinearity and stochasticity of the underlying biasing 
mechanism (Dekel, Lahav 1999). The deterministic linear bias <5 obj = & obj <5 maBB , for instance, implies that r obj = 1. 

Since the two-point correlation functions are not positive definite, however, the physical interpretations of the 
above definitions are somewhat ambiguous in that case. We will discuss this point in detail in Sec. 3. 4 where we 
relate them to those defined by the one-point statistics, i.e, & ob ] and r obJ -. In any case, this does not happen in our 
analysis below on scales of cosmological interest, and we adopt the above "conventional" definitions. In what follows, 
we focus on the scale-dependence, time-evolution of the above parameters, as well as on their behavior in redshift 
space. Before presenting the simulation results, we briefly summarize the existing theoretical models. They turn out 
to be useful in comparing and understanding our results despite that they are valid mainly in linear regimes. 



3.1. Theoretical models for biasing 
3.1.1. Peak model 

Density peaks are the first physical model of cosmological biasing proposed by Kaiser (1984), and then extensively 
studied by Bardeen et al. (1986). According to the latter, we consider the density peaks in the primordial random 
Gaussian field smoothed with a Gaussian window function. Then the height of each peak is defined as the ratio of 
the density contrast at the position smoothed with Rf and its RMS value, v = 5/ao. Although 5 and o~q in this 
context should refer to the values at some initial epoch, we take another (but equivalent) view that their values are 
linearly extrapolated to the present time. To be more specific, <jq is computed by setting I = in the following 
general expression for the Z-th order moment: 



of(Rf) 



^fc 2 <P(fc;^0)exp(-fc 2 i? 2 ), 



(9) 



where P{k; z = 0) is the linearly extrapolated power spectrum. 

Bardeen et al. (1986) showed that the number density of peaks, n poak , is given by 



1 



n poak (f) = 
where 

R* = V3 

7 = 
G{i,y) = 



" 2/2 G(7,7^), 



<Ji{Ri) 

<J2{RfY 
<7oORf)<720Rf)' 

exp[-(.T-y) 2 /2(l- 7 2 ) 



/>OC 

/ dx F(x) ■ 
Jo 



[2tt(1 - 7 2 )]!/2 
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erf 
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e -5x 2 /S + ' x 



-5x 2 /2 



(10) 

(11) 

(12) 
(13) 
(14) 



Mo, Jing & White (1997) derived a formula which relates the number density fluctuations of peaks to the mass 
density fluctuations <5 maS s at z. Denoting D(z) as the linear growth factor normalized to unity at z = 0, the 
extrapolated amplitude of mass fluctuations at the present time is S mass /D(z). Then, the fluctuation <5 poak with 
v > S^ B J[D(z)a ] is given by 



<U k (,, z) = (1 + S^) n ^ )l — ^ZIM _ ! ! * . 

n peak (z/) L D(z)a dv 



{^^^l 108 ^^}^' (15) 



in the limit of linear theory. Thus the peak model reduces to the scale-independent deterministic linear biasing: 



7 / \ "peak 

Opeak(^z) = 

"mass 



1 v 2 + .91 



9i 



D(z) a {R { y 
-iv 9G(7,y) 



G(7,7^) dy 



(16) 
(17) 
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Since our simulation data are for peaks with the peak height above v, we use the following effective biasing factor 
in the comparison below: 



- / p oak, off 



(> v,z) 



& P cak(V, z)n paak (is')dv' 



I 



n pee . k (i/)dv' 



(18) 



3.1.2. Halo biasing (Merging model) 

Another popular model of cosmological biasing is that one dark matter halo accommodate a single luminous 
object. This is established as a fairly realistic model for clusters of galaxies (e.g, Kitayama, Suto 1996), and may 
provide a reasonable approximation even for bright galaxies (Somerville et al. 2000). The number density of halos 
with mass M at z is well approximated by the Press-Schechter mass function (Press, Schechter 1974): 



n h&io (M,z)dM = 



1 P 



■ exp 



f^M A 3 {M,z) 
A 2 (M,z) = D\z) J 0fP(k-,z = 



A 2 (M,z) 

». 9 



dA 2 (M, z) 



dM 



dM, 



{kR^ 



[sin(fci? TH ) — kR TH cos(fci? TH )] 2 , 



(19) 
(20) 



where p is the mean mass density, S c is the critical density threshold of the spherical collapse, and the spherical 
top-hat radius R TH is (3M/47rp) 1/3 . 

An analytic model for dark halo biasing is constructed in a similar way as the peak model (Mo, White 1996; Mo 
et al. 1997), yielding 



<5 halo (M,z) = (l + 5„ 



^ n halo (M, z)\s c ->s c -s„ 
n halo (M, z) 



- 1 ~ jl - J- log [n halo (M, z)]| 6 n 



(21) 



Again in the above limit of linear regime, the halo biasing becomes linear, deterministic and scale- independent. In 
what follows, we adopt the fitting formula of Jing (1998) for & halo (M, z) which incorporates a correction factor for 
low mass halos: 



6 halo (M,z) = l^r + i 



0.5 



(0.06-0.02n eff ) 



1 



Z/h = 



n e s = 



A(M, z) 
d\nP{k) 



dink 



(22) 
(23) 
(24) 



k=2ir/R 

and introduce the effective halo biasing factor to compare with the simulation results: 



&halo,eff(> M,z) = 



/ 6 halo (M',z)n halo (M',z)dM' 

J M 

/•oo 

/ n halo (M',z)dM' 

JM 



(25) 



3.1.3. Time- evolution of biasing parameters in a number- conserving model 

As noted above, both peak and halo models result in linear, deterministic and scale-independent biasing on large 
scales. In this limit, time-evolution of the biasing parameters can be generally computed using the continuity 
equation if the number of objects is constant (number conserving model). This is indeed the case for the peak 
model, while not for the halo model. 

The explicit expression for the time-dependent biasing factor is obtained first by Fry (1996) in a context of the 
deterministic biasing model. This is later extended to the stochastic biasing case in the linear regime by Tegmark 
& Peebles (1998), and then in the weakly nonlinear regime by Taruya, Koyama & Soda (1999) and Taruya (2000). 
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We use the expression of Tegmark & Peebles (1998) in considering the evolution of peak bias, which relates (6 ob j, 
r obj ) at a given redshift z\ to those at an earlier epoch, z 2 (> z\): 

. , s V D 2 {*2, zi) - 2D(z 2 )b oh] (z 2 )r oh] (z 2 )D(z 2 , Zl ) + D 2 {z 2 )b ab3 {z 2 f 

D(zi) 

I \ _ -P(z2)frobj02)r obj 2 ) ~ gO^fl) . 07 n 
D{z 1 )b obl {z 1 ) 

where we define 



D(z 2 ,z 1 ) = D(z 2 )-D(z 1 ). (28) 
We note that in the case of a deterministic linear bias, r obj becomes unity and the above expression reduces to 

b obi ( Zl ) = l + ^[b obi (z 2 )-l}, (29) 

which is identical to the behavior of 6 pcak (eq.Jlq]) in a linear regime. 



3.2. Scale- dependence 

With the above theoretical predictions as a reference, let us examine the scale-dependence of the biasing param- 
eters. Figures 4 and 5 display the results for peak and halo models, respectively. Apart from the object-dependent 
amplitude of the parameters, the general trend is quite similar; on large scales, r obi is close to unity and the de- 
terministic linear bias is a reasonable approximation in this linear regime. In fact, both b p ^ k and b„2 approach 
the predictions ([l8|) and (|2^), respectively, as the separation r increases. On smaller scales, the biasing becomes 
nonlinear and stochastic (i.e., < 1) and the clustering amplitude relative to mass becomes generally stronger 

(bobj increases). Altho ugh bhaio levels off and then decreases below a few Mpc, this should be mainly ascribed to the 
exclusion effect due to the finite halo size (Mo, Jing, White 1997). 

Turn next to the redshift-space distortion effect on the scale-dependence. On large scales, the simulation results 
are understood fairly reasonably with a straightforward extension of a linear distortion model (Kaiser 1987) to 
the stochastic biasing model, which is described in Appendix; we obtain b„f^ b ab } and — ribj • On small 

scales, however, the degree of the redshift-space distortion depends sensitively on the nature of biased objects. Note 

(s) 

that the redshift-space distortion even leads to r nalo > 1, but this is simply due to the improper definition on small 
scales (see discussion in Section 3.4). We also find that while b p f ak is significantly suppressed with respect to b pc J k , 
the amplitude of b[f l0 is comparable to that of b na ^,. For the latter aspect, this implies that the pair- wise velocity 
dispersion, (v1 2 ob ,), on small scales is dependent on objects, and that (v{ 2 halo ) £ {v\ 2 maBB ) ~ {v\ 2 p Cak )- We directly 
compute the velocity dispersion and indeed confirm that the above relation holds in our simulation catalogues of 
peaks and halos. 



3.3. Time evolution 

Now we examine the time evolution of the biasing parameters. For this purpose, we consider the LCDM model 
and plot the biasing parameters as a function of the redshift in Figure 6. In order to probe the linear, quasi-nonlinear 
and nonlinear scales of the gravitational clustering, we select three separation lengths; r = 1, 5 and 23 /i -1 Mpc for 
peaks, and r = 7, 13 and 23/i~ 1 Mpc for halos. Again the qualitative behavior is quite similar in the peak and halo 
models regardless of the scale. Toward lower redshifts, r obj increases and approaches unity while 6 obj decreases. We 
plot the number-conserving model predictions (Tegmark, Peebles 1998) in Figure 6a using the values of simulations 
at z 2 — 3.4 for & poak (z2) and r pcak (z2) in equations j26[ ) and d27]). It turns out that the time evolution of biasing in 
peak model is well described by the number-conserving model even on fairly nonlinear scales. 

Figure 6b plots the results in halo model. Since our halo catalogues are selected by the mass, the number density 
of the selected halos changes with redshift according to the Press-Schechter mass function. Therefore the number- 
conserving model predictions do not match the results. Rather we plot the bias formula by Jing (1998) which seem 
to be in reasonable agreement with the evolution in simulations. Incidentally we create a different halo catalogue 
which conserves the number density by appropriately changing the selection mass threshold with z, and make sure 
that the evolution in such halo catalogues agrees well with the number-conserving model prediction. 
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3.4- Comparison with biasing parameters in terms of the one-point statistics 

As noted earlier, the biasing parameters defined in equations (fjj) and (||) are somewhat ambiguous on some scales. 
This leads to unphysical results on small scales that the cross-correlation factor r obj does not lie in the range [—1,1] 
and that the biasing factor 6 obj may become imaginary. Indeed, this behavior r abj > 1 shows up in the halo model 
(Figure 5). In the light of this, it is useful and instructive to relate our biasing parameters to those defined in terms 
of one-point statistics which are wel-defined on all scales. 

In general, the two-point statistics does not simply reduce to the one-point statistics because of the difference of 
their smoothing function. Nevertheless in the case of the top-hat window function, the variances of mass, objects 
and their cross-correlation can be explicitly written in terms of the corresponding two-point correlation functions as 
follows: 

/>2i?TH 

4(i? TH ) = 4tt / dr r 2 fo(r) F(r; R TH ), (30) 
Jo 

where the subscripts i and j mean either o or m. The filter function F(r;R TH ) is given by 
m R v 3 (r-2R TH ) 2 (r + AR TH ) 

F ^ R ^ = ^ (2R^f ' (31) 

and satisfies the normalization condition : 

/•2-Rth 

4vr / dr r 2 F(r;R TH ) = 1. (32) 
Jo 

Using the above relation, one can in principle reproduce the biasing and cross-correlation factors for one-point 
statistics in a straightforward manner. In practice, however, one has to take account of the finite resolution of the 
numerical data carefully. In particle simulations, the smallest scales are supposed to be dominated by the discreteness 
effect and the two-point correlation functions at r — are formally described by the sum of the Dirac delta function. 
In fact, this sensitively changes the amplitude of <Ty (_R TH ) on small scales. This implies that one has to introduce 
a particular model of two-point correlation functions at r = 0. For this purpose, we adopt the following procedure. 
We assume that the number density of each object on smallest scales is simply described by the sum of (randomly 
distributed) Dirac's delta function n obi (x) — ^2 u Sd(x — x n ). In this case, the Poisson process yields the following 
offset to the intrinsic clustering term cTj 2 i (i? TH ): 

^.trueC-ft™) = of 4 (i? TH ) + a 2 hot (R TH ), (33) 

where the first term in the right-hand-side is calculated from equation ( p0| ) using the numerical data of correlation 
functions above the scales of the resolution. The second term is given by 

-l 



a shot(^TH) — 



3 TI 



(34) 



in the case of the top-hat window function with n abJ being the mean number density of objects. 



Adopting equation (|33|), we estimate the variances untrue and then compute the biasing parameters b ob ] and r^j: 
bil] = (35) 



o- 2 



r£] = (36) 



> oo^ mm 



where a om = (<5 obj i5 maas ) is the one-point cross-correlation. 

In Figures 7 and 8, the biasing parameters for the one-point statistics are shown as a function of smoothing radius 
i? TH - On large scales, their behavior is similar to that for the two-point statistics, and the theoretical prediction 
is good agreement with the numerical results. On the other hand, on small scales, the point process dominates 
the clustering signal and scale-dependence of the biasing and cross-correlation factor is quite different from that of 
b obj and r abJ for both peaks and halos; ri b ] substantially deviates from unity and & ob ] becomes quite large as the 
smoothing radius R TH decreases. The effect becomes more significant in the redshift-space and/or at higher redshift, 
since the redshift-space distortion and redshift evolution tend to weaken the intrinsic clustering amplitude. 

The comparison of those biasing parameters between the one- and two-point statistics suggests that our results 
are reliable on scales larger than ~ lft, _1 Mpc for peaks and ~ 5/i -1 Mpc for halos where the small-scale modeling 
does not seriously affect the biasing parameters. 
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4. Conclusions and discussion 

We have extensively examined the scale-dependence, time-evolution, and redshift-space distortion of biasing in 
peak and halo models, with particular attention to the nonlinear stochastic nature. Especially, we quantify the 
biasing properties defined in terms of the two-point statistics and the results are compared with those in the one- 
point statistics. Our main results are summarized as follows: 

(1) We have quantitatively demonstrated the degree of nonlinearity and stochasticity of biasing in peak and halo 
models. In particular, biasing is significantly nonlinear and stochastic in nonlinear regimes of gravitational clustering, 
due to the gravitational evolution for peaks and volume exclusion for halos, which results in strong scale-dependence 
of the biasing parameters on small scales ( < a few Mpc) . 

(2) On scales larger than 10/i _1 Mpc, however, a linear deterministic and scale-independent biasing remains a 
reasonable approximation in peak and halo models. While the biasing parameters in this linear regime evolve as a 
redshift, their evolution is consistent with the existing theoretical model predictions. 

(3) Redshift-space distortion of the biasing on small scales is quite different in peak and halo models because the 
pair-wise velocity dispersions are quite sensitive to those objects. This difference affects the stochasticity of biasing 
in redshift space. The linear redshift-space distortion, on the other hand, is very similar and explained by extending 
the model of Kaiser (1987) to the stochastic biasing case (see Appendix). 

The biasing models considered here are rather idealistic for the observed luminous objects which form through the 
complicated processes including the gas dynamics and the radiative transfer. Nevertheless the peak and halo models 
contain the most important aspect of the formation of astrophysical objects, the nonlinear gravitational growth and 
the subsequent merging, in the scenario of hierarchical galaxy formation. This is why our results are qualitatively 
similar to the recent numerical work of Colin et al. (1999) and Somerville et al. (1999), and also to predictions of 
an analytical model of Taruya & Suto (2000). Moreover cosmological simulations of galaxy formation by Yoshikawa 
et al. (2000) indicate that the one-to-one correspondence between halos and galaxies may be a reasonably good 
approximation on large scales and the nonlinear clustering of peaks traces the galaxy distribution fairly well even 
on small scales. 

The validity of linear deterministic and scale-independent relation on large scales may be a natural consequence of 
the hierarchical clustering scenario (e.g., Matsubara 1999). Such scale-independence probably comes from the fact 
that the statistical feature of large-scale clustering of peaks and halos is simply characterized by a single parameter; 
the mass threshold M for halos and the density height v for peaks. Our important result is that the relation 
S obi ~ b obj d maaB still holds in quasi-linear regime. If this is the case, the interpretation of clustering of realistic objects 
is straightforward. On the other hand, the scale-dependence of the biasing becomes crucial in the nonlinear regime 
of density fluctuations. These should be kept in mind in interpreting the future redshift surveys of galaxies. 

We thank Ue-Li Pen, for pointing out the importance of the relation between the one- and two-point statistics 
(eq.p(J]), which improves our analysis and leads to the discussion in Section 3.4. Y. P. J. and A.T. gratefully 
acknowledge support from a JSPS (Japan Society for the Promotion of Science) fellowship. Numerical computations 
were carried out on VPP300/16R and VX/4R at ADAC (the Astronomical Data Analysis Center) of the National 
Astronomical Observatory, Japan, as well as at RESCEU (Research Center for the Early Universe, University of 
Tokyo) and KEK (High Energy Accelerator Research Organization, Japan). This research was supported in part 
by the Grant-in- Aid by the Ministry of Education, Science, Sports and Culture of Japan (07CE2002, 12640231) to 
RESCEU, and by the Supercomputer Project (No.99-52, No.00-63) of KEK. 



Appendix. Linear stochastic biasing in redshift-space 

This appendix describes the generalization of the the linear redshift-space distortion model of Kaiser (1987) in 
the context of the stochastic biasing. Although we consider the biasing for the one-point statistics, the results are 
also applicable to the biasing for two-point statistics as long as the correlation functions are positive. 

First recall that the density contrast of dark matter and objects in k-space is given by 




ei(k)+/M 2 ^i(k), 
^(k)+v4i(k), 



(Al) 
(A2) 
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where [i = k\Jk is the direction cosine in k-space (fcy is the comoving wave number parallel to the line of sight), and 



/ 

Q(z) = 
X(z) = 



d1nD(z) 
din a 



Sl(z) 



0.6 



A(z) 
70 



1 



Q(z) 



Hp 
H(z)\ 

Ho 



H(z) 

Equations flAl| ) and (f 

(€L€1)(m) = 



A . 



(A3) 
(A4) 
(A5) 



relate the auto- and cross-correlation functions in redshift space to those in real space: 
(l + 2/ / , 2 + /V 4 )(eL^i>, (A6) 



1 



(R)x(R)| 



.,(R)JR) 



+ /M 2 1 + b[ 



>,(RUR) 



(A7) 
(A8) 



where the brackets indicate the ensemble average, (3 
and cross-correlation factors in real space: 



f/b^j and the quantities and r^i denote the biasing 



(R) 



„(R) - 



(R) r(R) 



\<«2 



(A9) 



?(R)\ 



(A10) 

The biasing and cross-correlation factors in redshift space are defined similarly using equations (A6)-(IaE). After 



cross-con 
integrating over /i, we find 

b (S) 



l(R) 



\ 



„(S) _ 



(i + f/frg + 1^) 
(i + !/ + |/ 2 ) ' 

rS ) + ^(l + 6SVS , ) + */i9 



bj ^o bj 



(! + !/' 



In the deterministic case r, 



^/ 2 )(1 
1, one can show that 



¥ 2 ) 



dS) 

J obj 

AS) 



< 
< 



■>obj 1 



1, 



(All) 



(A12) 



(A13) 
(A14) 



from equations (All) and (A12). The inequality (A13) is understood by the squashing effect which is more appre- 
ciable for 6(f2>B if^Ibj > <^SL (see eqs.[Al] and A2|). The inequality ( A14| ) merely reflects the fact that the random 
peculiar velocity generates an additional scatter around the mean biasing relation in redshift space. 

( R) A k 

In a more general case, the biasing is not deterministic i.e., r ohi < 1. Even then the inequality (A13) still holds as 



long as bi^ > 1 and r ab ]' < 1. The inequality (A14), on the other hand, no longer holds, and in fact, the opposite 
inequality seems to be a fairly general result; Figure 9 plots r„ h ] — as a function of r^j for different values of 
/; / = 0.3 (solid lines), / = 0.6 (dashed lines), and 0.9 (dotted lines) with b obj = 3.0 (left-panel) and b obi ' = 5.0 



„(«) 



(right-panel). Except for a small region around r, 



1. 



AS) 



is generally larger than 



„(«) 



It seems that the 



stochasticity is effectively reduced due to the redshift-space distortion as increasing the stochasticity in real space. 
Moreover this tendency becomes stronger when the factor / becomes large, or equivalently, the density parameter 
Qo approaches unity. 



( S) ( R) 

Thus in linear theory of redshift-space distortion, we conclude fairly generally that b ab - < 6, ' 



long as a certain degree of the stochascity exists in real space (i.e. 



bj and r obj 



(R) 
K hj as 



< 1) as Figure 9 exhibits. 
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Table 1. Simulation model parameters. 



Model 


n 


A 


rt 


CT8 


N 




realizations 


SCDM (Standard CDM) 


1.0 


0.0 


0.5 


0.6 


256 3 


4.5 x 10 11 


3 


LCDM (Lambda CDM) 


0.3 


0.7 


0.21 


1.0 


256 3 


1.3 x 10 11 


3 


OCDM (Open CDM) 


0.3 


0.0 


0.25 


1.0 


256 3 


1.3 x 10 11 


3 



t the shape parameter of the power spectrum. 
* mass of a single dark matter particle. 



Table 2. The peak catalogue. 





SCDM 

Art 


b* 


LCDM 

N 


^pcak,cff 


OCDM 

N 


^pcak,cff 


1.0 


6.2 x 10 5 


1.2 (1.8) 


5.3 x 10 5 


1.2 (1.5) 


5.5 x 10 5 


1.2 (1.4) 


2.0 


2.6 x 10 5 


1.6 (2.8) 


2.2 x 10 5 


1.5 (2.2) 


2.3 x 10 5 


1.5 (1.9) 


3.0 


5.0 x 10 4 


2.0 (4.1) 


4.5 x 10 4 


1.8 (3.1) 


4.6 x 10 4 


1.8 (2.6) 



t average total number of peaks among three realizations. 
* effective biasing factor of peaks at z = 0.0 (z = 2.2). 



Table 3. The halo catalogue. 



Mth^Mo] 


SCDM 

JVt 


"halo, oft 


LCDM 

N 


&halo,eff 


OCDM 

N 


&halo,eff 


2 x 10 12 
5 x 10 12 
1 x 10 13 


- (- -) 
1.2 x 10 5 (1.4 x 10 4 ) 
5.7 x 10 4 (2.3 x 10 3 ) 


1.1 (5.2) 
1.1 (5.2) 
1.3 (6.7) 


7.3 x 10 4 (3.9 x 10 4 ) 
3.0 x 10 4 (1.1 x 10 4 ) 
1.5 x 10 4 (3.3 x 10 3 ) 


0.9 (2.6) 
1.0 (3.3) 
1.2 (4.1) 


8.2 x 10 4 (6.6 x 10 4 ) 
3.4 x 10 4 (2.1 x 10 4 ) 
1.7 x 10 4 (8.3 x 10 4 ) 


0.9 (1.8) 
1.0 (2.3) 
1.2 (2.8) 



t average total number of halos among three realizations at z = 0.0 (z = 2.2). 
* effective biasing factor of halos atz = 0.0(z = 2.2). 
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(a) Real space Redshift space z = o.O 




— 300 h^Mpc — 



(b) Real space Redshift space z = s.S 




— 300 h^Mpc — 



Fig. 1.. Top-view of distribution of objects in real (left panels) and redshift (right panels) spaces around the fiducial observer at the 
center; dark matter particles (top panels), peaks with v > 2 (middle panels) and halos with M > 1.3 X 10 12 Mq (bottom panels) in 
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Fig. 2.. Auto- and cross-correlation functions of dark matter and peaks in SCDM (left panels), LCDM (middle panels) and OCDM 
(right panels). Different symbols indicate the results in real space (open squares for v > 3, filled triangles for v > 2, open circles for 
v > 1, and crosses for dark matter), while different curves indicate those in redshift space (dashed for v > 3, dot— dashed for v > 2, 
solid for v > 1, and dotted for dark matter), (a) z = 0, (b) z = 2.2. 




Fig. 3.. Same as Fig. 2 but for halo model; open squares and dashed lines for M > 10 13 h~ 1 MQ, filled triangles and dot-dashed lines 
for M > 5 X lO 12 fc^" 1 M0, open circles and solid lines for M > 2 X 10 12 /i — 1 Mq, and crosses and dotted lines for dark matter. For 
SCDM model, we only plot the correlation functions with M th = 5 X 10 12 , 10 13 h _1 M Q . (a) z = 0, (b) z = 2.2. 
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Fig. 4.. Scale-dependence of the biasing and cross-correlation factors for two-point statistics of peak model. Different symbols and lines 
have the same meaning as in Fig. 2. The arrows in the bottom panels indicate the model predictions of Mo, Jing, & White (1997) in 
real (dotted) and redshift (solid) spaces, (a) z = 0, (b) z = 2.2. 
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Fig. 5.. Same as Fig. 4, but for halo model. Different symbols and lines have the same meaning as in Fig. 3. The arrows in the bottom 
panels indicate the model predictions of Jing (1998) in real (dotted) and redshift (solid) spaces, (a) z = 0, (b) z = 2.2. 
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Fig. 6.. Time evolution of the biasing and cross-correlation factors in LCDM model, (a) Open (filled) symbols indicate the peak 
simulation results for v > 3 (squares), v > 2 (triangles), and v > 1 (circles) peaks in real (redshift) spaces. Dotted and solid lines 
correspond to the model predictions of Tegmark & Peebles (1998) in real and redshift spaces, respectively, normalized by the data 
at z = 3.4. (b) Open (filled) symbols indicate the halo simulation results for M > 10 13 A _1 M Q (squares), M > 5 X 10 12 /i _1 M Q 
(triangles), and M > 2 X 10 12 /i _1 Mq (circles) peaks in real (redshift) spaces. Dotted and solid curves correspond to the model 
prediction of Jing (1998) in real and redshift spaces, respectively, (b) halo model. 
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Fig. 7.. Scale-dependence of the biasing and cross-correlation factors for one-point statistics of peak model. Different symbols and lines 
have the same meaning as in Fig.2. The arrows in the bottom panels indicate the same model predictions as in Fig. 4. (a) z = 0, (b) 
z = 2.2. 




Fig. 8.. Same as Fig. 7, but for halo model. Different symbols and lines have the same meaning as in Fig. 3. The arrows in the bottom 
panels indicate the same model predictions as in Fig.5. (a) z = 0, (b) z = 2.2. 
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Fig. 9.. The difference r^J - r^J as a function of r$ for / = 0.3 (solid), / = 0.6 (dashed), and 0.9 (dotted) with b ( a 
(left-panel) and b^j =5.0 (right-panel). 



